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We introduce and analyze a lattice model of anyons in a periodic potential and 
an external magnetic field which exhibits a transition from a Mott insulator to a 



^ ■ quantum Hall fluid. The transition is characterized by the anyon statistics, a, which 



can vary between Fermions, a = 0, and Bosons, a = 1. For bosons the transition 
is in the universality class of the classical three-dimensional XY model. Near the 
Fermion limit, the transition is described by a massless 2 + 1 Dirac theory coupled to 
a Chern-Simons gauge field. Analytic calculations perturbative in a, and also a large 
N-expansion, show that due to gauge fluctuations, the critical properties of the tran- 
sition are dependent on the anyon statistics. Comparison with previous calcualations 
at and near the Boson limit, strongly suggest that our lattice model exhibits a fixed 
line of critical points, with universal critical properties which vary continuosly and 
monotonically as one passes from Fermions to Bosons. Possible relevance to experi- 



ments on the transitions between plateaus in the fractional quantum Hall effect and 
the magnetic field-tuned superconductor-insulator transition are briefly discussed. 
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I. INTRODUCTION 

A very powerful approach for treating strongly correlated quantum models in two- 
dimensions has been to transform the statistics of the particles by attaching statistical flux 
tubes |]]]. In this way, a fermion model can formally be transformed into a boson model and 
vice versa. Moreover, models describing particles with fractional statistics (anyons) can be 
mapped either way, into a fermion or a boson model. This mapping is usually accompanied 
by a flux smearing 'mean field' treatment, originally introduced by Fetter, Laughlin and 
Hanna ||. In this mean field approach the statistical flux tubes are 'detached' from the 
particles and smeared uniformly in space. Taken together with the above mapping, this 
mean field approach effectively 'trades in' statistics for an external magnetic field! This 
approach serves as the key ingredient in various recent theories of the quantum Hall effect: 
For example the so-called Ginzburg-Landau approaches 0-fl which focus on an underlying 
Bose condensation, the hierarchical construction due to Jain |||7| which relates integer and 
fractional Hall states, and most recently to theories of the half-filled Landau level || as a 
pseudo-fermi liquid. It has also recently been applied in attempt [[J to relate the continuous 
transitions between plateaus in the fractional Hall effect to the transition between plateaus 
in the integer effect. 

It is generally believed that this flux-smearing mean field approximation should be legit- 
imate when the mean field state is an incompressible fluid with a gap, such as a quantum 
Hall fluid. Indeed, various numerical exact diagonalizations |I| on small systems appear 
to confirm this belief. However, when the reference state is gapless, as for example at a 
transition between plateaus in the quantum Hall effect or in the fermi-liquid theories of the 
half-filled Landau level, the approximation is much more highly suspect. In particular vari- 
ous conclusions arrived at in Ref. || which relate detailed critical properties at transitions 
between different quantum Hall plateaus are worth detailed scrutiny. 

The purpose of the present paper is to introduce the simplest possible anyon model 
which exhibits a continuous phase transition, and analyze in detail the universal critical 
properties near and at the transition. The model is characterized by the statistics of the 
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particles, which can be tuned continuously from Bose to Fermi. The particles are subject to 
a periodic potential with a period commensurate with the particle's density, corresponding 
to one particle per unit cell. Moreover, the model includes an external magnetic field, 
taken as zero in the boson limit, one flux quantum per particle in the fermion limit, and 
proportional to the statistics in the anyon case. In this way a 'mean-field' treatment of 
the model, as described above, would result in behavior fully independent of the particle's 
statistics. Controlled calculations of the actual critical properties then allow for a direct 
check on the validity of the mean-field approach. 

In the boson limit, the model exhibits a Mott-insulator to superfluid transition which is in 



the universality class of the classical 3D XY model [|ll|] . Critical properties can be extracted 



via Monte Carlo simulations or 1/N expansions |TI|. As we show below, in the Fermion 
limit the model is soluble and exhibits a 'gap-closing' transition between a band-insulator 
and an integer quantum Hall state. We analyze the critical properties in this case, which are 
described by a massless 2+1 Dirac equation, and find that they are most certainly different 
from the 3D XY model, which reveals the inadequacies of the 'mean-field' approximation. 

For anyon statistics the model exhibits a transition between a Mott-insulator and a quan- 
tum Hall fluid. Can the critical properties be extracted for this anyonic Mott transition? 



As shown recently by Wen and Wu [|13l , it is possible to perform a renormalization group 
(RG) calculation perturbative in deviations of the particle's statistics from the bosonic point. 
Within a controlled 1/N expansion they find that the critical properties vary continuously 
with the particle statistics! In RG terminology they have found a 'fixed line' parameterized 
by the particle's statistics. In this paper we perform a similar perturbative analysis ex- 
panding around the Fermi end, described by the Dirac equation. We likewise find that the 
critical exponents and other universal properties vary continuously with statistics. Again, 
this statistics-dependence is shown to originate from gauge-field fluctuations near the criti- 
cal point. Moreover the sign of our perturbative results for the critical exponents near the 
Fermion point, suggest that the exponents vary monotonically upon moving along the fixed 
line between the Fermion and Boson transitions. We suspect that this might be a generic 



feature of quantum phase transitions involving fractional statistics particles. This would 
imply that the leading perturbative corrections to the flux-smearing mean field treatment 
approach monotonically the exact behavior. 

The dependence of the detailed critical properties on the statistics of the underlying 
particles in our model system, strongly suggests that, in contrast to recent speculations 
||, the phase transition between integer plateaus in the quantum Hall effect is not in the 
same universality class as that between fractional plateaus. Indeed, in view of our results, 
it would be much more natural to argue that the exponents will depend on the statistics of 
the condensing quasiparticles, which of course depend on which two fractions the transition 
is between. However since both disorder and long-ranged Coulomb interactions are relevant 
perturbations at the Mott transition studied here, we cannot make any definitive statements 
about the universality class appropriate to experimental quantum Hall phase transitions. 

The paper is organized as follows. In Section II we introduce the lattice anyon model, 
and discuss the continuum limit, with special emphasis on the fermion limit. The resulting 
critical theory is a massless 2+1 Dirac equation coupled to a Chern-Simons term. Section 
III is devoted to a detailed RG analysis of this fermionic Chern-Simons theory, perturbative 
in deviations of the statistics away from Fermi, as well as in a large N limit. In addition 
to the critical exponents v and 77 we calculate the universal Hall and longitudinal conduc- 
tivities. Finally, Section IV is a discussion section with emphasis on possible relevance to 
experimental systems. 

II. THE MODEL 

In this Section we introduce our lattice anyon model, and in the fermion case extract the 
appropriate continuum limit. 

The Hamiltonian we study is defined for convenience in terms of spinless fermion oper- 
ators: 




(1) 



where the fermion operators satisfy the usual anticommutation relations: 
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{ Cl ,c]} = S t3 . (2) 

Here i and j label sites of a square lattice, tij is a hopping matrix element and Sj = +1 for 
sites of one sub-lattice and —1 on the other sub-lattice. The second term in the Hamiltonian 
is thus a staggered potential with strength m. We restrict attention exclusively to half-filling, 
so that the number of particles equals the number of sites on either of the two sub-lattices. 

An applied magnetic field, B, and flux tubes attached to the fermions both enter via a 
gauge field: ty = \tij\e lAiJ . For later convenience we take the hopping strength \tij\ to be 
equal to t for nearest neighbor sites, t' /4 for next nearest neighbors and zero otherwise. The 
gauge field is chosen so that 

V x A = B + ap , (3) 

where Vxi denotes a lattice curl, that is an oriented sum of around plaquettes, and B is 
a uniform external field in units of a flux quantum, h/e, per elementary square of the lattice. 
Here p, which is an operator, denotes the particle density so that Eq.(3) must be thus taken 
as a constraint. The statistics parameter a gives the strength of the statistical flux tubes 
which are thereby attached to each particle. With a = 1 the particles are transmuted into 
bosons, whereas anyons correspond to a between zero and 1. Finally, to fully specify the 
model we choose the external field as B = (1 — a)/2, so that at half-filling in a flux-smearing 
mean-field treatment the Hamiltonian reduces to bosons in zero field, for all a. 

Notice that the model has been constructed so that in the boson limit (a = 1) it is 
formally equivalent to a lattice model of hard core bosons in zero magnetic field. Moreover, 
at half-filling the boson density is commensurate with the period of the staggered potential. 
This boson model is expected to undergo a Mott-insulator to superfluid transition Jll] 



as 



the ratio of the kinetic energy (t, t') to the staggered potential (m) is varied. This zero 
temperature quantum transition is expected to be in the universality class of the classical 
3-dimensional XY model. The appropriate coarse-grained continuum theory to describe the 
behavior near the transition is then simply a <ft A complex scalar field theory. 
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For a near 1, the model corresponds to anyons with statistics 'close' to bosonic, in a 
weak magnetic field. In this case the model is expected to exhibit a transition from a 
Mott-insulator to a gaped quantum Hall fluid. The transition can be studied by adding 
a Chern-Simons term to the scalar field theory, which attaches (1 — a) flux tubes to each 
boson. The corresponding critical behavior has been studied in a recent large- N approach 
13| , and the exponents found to vary continuously with a. 

When a = the Hamiltonian is simply that of non-interacting spinless fermions in a 
magnetic field, and can be solved by straightforward diagonalization. Specifically, with 1/2 
flux per square and the staggered potential the unit cell has four sites fllijl , as sketched in 



Fig. 1. Extracting the band structure thus involves diagonalizing a 4 by 4 matrix. Doing so 
reveals that generally there are two bands, symmetric about zero energy (half-filling) , with 
a band gap about the origin. Depending on the relative sizes of t' and m the bands can 
be either 'Landau levels', which contribute a Hall conductance, or conventional insulating 
band with no Hall effect [[HJ. Of interest to us is the transition between these two phases, 
where the band gap closes. (Across the transition the Hall conductance jumps from unity to 
zero, in units of e 2 /h with e the electric charge of the fermion.) The band structure reveals 
that the gap closes at a single point in crystal momentum (k-) space. By focusing on values 
of k near this point, we now construct a continuum field theory for the behavior near the 
transition. 

To this end, we choose an explicit gauge for the A^-, which is shown in Figure (1). 
Moreover, we put t—1 and take t' and m much smaller than 1. In this limit the band gap 
vanishes when t' = m and occurs at (it, tt) in the Brouillon zone. Upon linearizing the 
momentum about this point the Hamiltonian can be cast in the form: 

H = [ £ ci(p)H ab (p)c b (p) , (4) 

J P 0.6=1 

where c a (p) denotes the Fourier transform of the electron operators on the 'a' sites (see 
Fig. 1) at momentum (n,7c) + p. The 4 by 4 matrix H ab can be written 



Hab{p) = {PyOy ~ m(T z ) <g> 1 + (p x (J x + t'(T z ) <g> T y , (5) 



where <r M and r M denote the usual Pauli matrices and <8> is a direct product. Upon performing 
a basis rotation on by 90 degrees about the x-axis, the Hamiltonian can be transformed 
into a 2 by 2 block diagonal form, 

H ab (p) = H + (l + r g )/2 + if_(l - r z )/2 , (6) 

with 

H± = ±p x cr x + py<T y + (±t' - m)a z , (7) 

Notice that H + is a two-dimensional Dirac Hamiltonian with mass, M — t' — m, which 
vanishes at the transition. Under a x — *> —a x , H_ is also a Dirac Hamiltonian except with 
mass t' + m. This mass remains non-zero when m = t', and thus non-critical. As far as 
critical properties near the transition are concerned we can safely ignore the massive field 
and focus exclusively on the single Dirac field which goes massless. Since this theory is 
non- interacting we can easily extract all of the relevant critical properties (see Section 3). 

When the statistics parameter a is non-zero, but small, progress can be made by ex- 
panding around the fermion point (ie. the Dirac equation). The appropriate continuum 
theory near the transition in this case can be simply obtained from the Dirac equation by 
minimal coupling to the gauge field whose associated flux must be attached to the fermions. 
This can be achieved in the usual way by the addition of a Chern-Simons term |l7j to the 
2+1 Dirac Lagrangian. In the next Section we perform a renormalization group analysis 
on this 2+1 Dirac plus Chern-Simons theory to extract critical properties near the anyon 
Mott transition. Two expansions are available for calculating them: the weak 'coupling' 
expansion, in which the small parameter is a, the deviation from Fermi statistics, and the 
1/N expansion, if one consider the case with N species of anyons. 

III. CRITICAL PROPERTIES 

The quantities of primary physical interest which characterize the Mott transition are 
the transverse and longitudinal conductivities at the critical point, and the usual critical 



exponents 77 and v. It has recently been argued [0] that in general at zero temperature 



quantum phase transitions in two spatial dimensions, the conductivity should be universal, 
and we indeed verify this below for the anyon Mott transition. As usual we define the 
exponent v in terms of the correlation length which diverges as M _I/ , upon approaching the 
transition by taking the Dirac mass M = m — if to zero. We define the exponent rj via the 
decay of the two-point anyonic correlation function 

< ^(x)i>(0) >~x~ {1+v) , (8) 

where ip is the anyon field operator in either a fermionized or bosonized representation. 
Implicit in this definition is the assumption that the anyons are converted into Fermions or 
Bosons by attaching flux tubes in the Landau gauge, where the vector potential is divergence 
free. Notice that in the above definition of r\ we have pulled out a factor of 1, so that it 



coincides with the usual exponent rj at the Boson Mott-insulator to superfluid transition |TT 
where the 1 is twice the canonical dimension of the boson field and ij twice the anomalous 
dimension. However, as we see below, at the Fermion Mott transition the scaling dimension 
of the Fermi field is 1, and there is no anomalous dimension, so that in this case rj — 1, 
rather than twice the anomalous dimension. 

To extract critical properties near the Fermion point we employ a renormalization group 
(RG) analysis in the framework of continuum Euclidean field theory by applying 

standard diagrammatic perturbative techniques [^1 . 

1) Field Theory Approach 

As shown in Section II, the long-length scale and low energy physics near the anyon Mott 
transition can be described by two D=2+l Dirac fermions, one of which becomes massless 
at the transition, which are coupled to a Chern-Simons gauge field a^. Specifically, from 
the Hamiltonian H in Eq. (|6|) and the constraint in Eq. ([3D, the appropriate Euclidean 
Lagrangian is simply 

L = X)V4b7i(d/j - W^) + iM±]i/;± + i^e^xa^ax ■ (9) 
± 1 

where ip± are the two species of Dirac fields. Here the Dirac matrices in three dimensions 
are 7 M = ia^, with (/i = 1, 2, 3) Pauli matrices, so that 
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rf Y = -5^1 - e^V , Tr(7^) = . (10) 

We have normalized the coefficient of the Chern-Simons term, and the 'coupling constant' 
g 2 is essentially the anyon statistics, a, measured from the fermion point: 

g 2 = 2na . (11) 

A source electromagnetic field which can be used to calculate the conductivites within linear 
response is minimally coupled to the Dirac fields in the usual way. It is worth emphasizing 
that the sign of the Chern-Simons term in the Lagrangian (9) must be chosen correctly 
in order that the flux tubes attached to the particles are in the opposite direction to the 
external physical magnetic field. 

The phase transition between the quantum Hall phase and the Mott insulator occurs 
when one of the Dirac masses, say M = M + , passes through zero. Before discussing the 
associated critical properties it is instructive to first recover the expected behavior of the 
two phases from the above Lagrangian. Provided we put in an ultraviolet cutoff, which is 
appropriate for the original lattice theory, when both masses are non-zero straightforward 
perturbation theory in powers of g 2 is convergent. Consider, for example, the fermion 
'polarization tensor', defined as the Fermion current- current correlation function, which we 
denote as Yl^ v {p). Current conservation (or gauge invariance) implies that the polarization 
tensor is always transverse: p^JX^ip) = 0. Such a tensor in D = 3 can be decomposed into 
an even and odd part: 

LV(p) = U e (p)p(S^ - ^) + U (p)e^ xPx . (12) 

where we have used the notation p = \p\. From standard linear response theory, the two 
parts are related to the Fermion longitudinal and transverse conductances, in units of e 2 /h 
with e the charge, by 

&f,xx = --II e (0) , a F ,xy = --n (o) . (13) 

a a 

The conductance of the anyons has an additional contribution coming from the attached 
flux tubes. Formally this is due to the difference between the anyon current operator, 
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^^(dfj, — iga p )ijj, which couples to the physical electromagnetic field, and the Fermion 
current operator: r y li ^)*d li ij}. (see Fig. 2). As can be easily shown (see e.g. ||13||), the anyon 
resistivty tensor, pij, is simply related to the fermionic resistivity tensor via: 

Pij = PF,ij ~ tijOi. (14) 

Previously there has been an intense effort in analyzing the Chern-Simons gauge theory 
coupled to a Dirac field in perturbation theory employing standard diagrammatic techniques 
j. Specifically, for a massive Dirac field it has been shown that the one- 



loop diagram in Fig. 3a contributes a finite contribution to the odd part of the Fermion 
polarization, H^(p = 0) = +(a/2)sgn(M±). Notice the dependence on the sign of the 
Dirac mass. Moreover it has been proved p3| , p6| that in the massive theory there are no 
additional non-zero corrections to the Fermion polarization, n e o (0), to all orders in a. Thus 
away from the Mott transition the Fermion conductivity tensor which follows from Eqn. |13| 
is given by 

= -tUj s 9n{M±). (15) 
1 ± 

In the Mott insulating phase, where the two masses have different signs, the conductivity 
tensor vanishes as expected, whereas in the Quantum Hall effect phase one has simply <jF,ij = 
—€ij. The total anyon resistivity tensor in the quantum Hall phase is then p^ = e^(l — a). 
Notice that when a = 1, which corresponds to the bosons superfluid phase (in zero magnetic 
field), the Hall resistivity vanishes as expected. 

We now turn to the critical properties near the transition where the mass of one Dirac 
field, M = M + vanishes. In the following we focus exclusively on the massless Dirac field, 
since the contributions from the massive field are known to all orders, as described above. 
Consider a renormalization group (RG) transformation for this massless field, which involves 
integrating over a shell of (three-)momenta p in a shell between A and A/6, with b > 1. 
To complete the RG we rescale the momenta by b and the fields ip and a by b 1+1 ^ and 
6 1+7a , respectively. The anomalous dimensions will be chosen to keep the coefficients of the 
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quadratic terms in the Lagrangian Eqn. (|9|) fixed. Power counting at the Gaussian (g = 0) 
fixed point reveals that the coupling constant g is indeed dimensionless. A straightforward 
perturbative renormalization group can then be carried out using standard methods. As 
defined above in Equation (|J) the critical exponents 77 and v are related to the anomalous 
dimension of the massless Fermion field, ip, and that of the composite operator ipip, which 
we denote by 7^ [fL9| , pO|l : 



77 = 1 + 27,/,, v 1 = l->y^ (16) 



It can be easily shown ( |26[) at one-loop order that the beta-function, (3(g) = 
—dg/dln(b) = 0(g 5 ), vanishes. Indeed it has been suggested |26j and shown explicitly 
to second order, that the beta-function vanishes identically: 

(3(g) =0. (17) 

This implies that there is a line of fixed points or critical points, parametrized by the anyon 
statistics a\ Unfortunately all of the anomalous dimensions vanish at one-loop order, so 
it is necessary to go to second order to evaluate the leading non-trivial corrections to the 
exponents. In addition to the critical exponents, the universal anyon conductivity at the 
critical point is of interest. As first emphasized by Semenoff, Sodano and Wu ||26|| , in contrast 



to the massive theory, in the massless theory there is a non- vanishing two loop contribution 
to the conductivity. Indeed, the exact value of this correction was obtained subsequently by 
Chen p9|j30[| . In the following we will give a unified description of the two-loop perturbative 
results for the critical properties of the model (|9|) in the context of the anyon Mott transition, 
together with some new calculations when the quantities of interest are not available in the 
literature. 

2) Weak- Coupling Expansion 

In this subsection we carry out a two-loop perturbation expansion in powers of the 
coupling constant g 2 to extract critical exponents and the universal conductivity. We take 
the Landau gauge, by adding a gauge fixing term, (l/2\)(d fl a^) 2 , to the Lagrangian @ 
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and then take the limit A — > in the resulting propagator for the gauge field. With the 
normalization chosen as shown in eq.(^) the appropriate Feynman rules for the massless 
Dirac field in the Landau gauge are 

the fermion propagator S (p) = — - , (18) 

i p 

the gauge propagator Gn U (p) = 7, — , (19) 

pi 

the interaction vertex ^oip) = Wl 1 * ■ (20) 

The advantage of the Landau gauge is that no infra-red divergences would appear pertur- 
batively in this gauge |2|j2|^8|. 

Since we require a two-loop calculation, rather than regularizing in the ultraviolet with 
a finite momentum cutoff, A « 1/a, it is more convenient to regularize by dimensional 
reduction and also employ a minimal subtraction method. In this approach, which was first 
suggested and extensively used in Refs. [^^HJ for the D = 3 Chern-Simons gauge theory, 



all momentum integrals in the loop-expansion are continued to general dimension D: 

,H /'•"' " /' 7^ ■ (21) 



(2tt) 3 ^ J (2tt) d 

where, in order to balance the dimension, one must introduce a parameter /1 which has 
dimension of mass. All vector, tensor and spinor quantities and in particular the symbol 
e iiu\ faofc appears in the original Feynman integrands, though, are always treated as if they 
were formally three dimensional. This implies that we will use the identities 

5 W = 3 , (22) 
e^V rr? = 5 UT 5 Xri - 5 uri 5 XT (23) 

before performing any momentum integrals. Finally, we adopt a minimal subtraction method 
for renormalization, by removing the simple poles in e = 3 — D |34|] and setting all higher- 
order poles to be zero. The parameter /1 we introduced in Eqn. ( pip represents, as usual, 
the "renormalization point" in the minimal subtraction approach. In all cases that have 



been checked |27H30l, including pure Chern-Simons (Abelian or non-Abelian) gauge theory, 
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coupled or not to massless or massive boson or Fermion matter fields, this approach to 
regularization and renormalization has been shown to respect both gauge invariance and 
Ward identities, at least up to two loops. 

We consider first the universal conductances at the critical point (M = 0), focussing 
exclusively on the contribution from the massless Dirac field. (The contribution from the 
massive field is given in Eqn. (^).) The one-loop polarization tensor IT^Qo), due to the 
(massless) fermion bubble (Fig. 3a), is purely symmetric in the indices /i and v and given 
by: 

n«(p) = -^, n«(p) = o. (24) 

Thus upon using Eqn. (|Tl~D and fll3|), to leading (zero'th) order in a, the total Fermion (and 
anyon) conductivity at the Mott transition is given by ap, X y — 1/2, coming solely from the 
massive Dirac field, and (Jf,xx = n/8. It is interesting that this value is precisely equal to the 
value of a xx obtained at the Boson Mott-insulator to superfluid transition within a large-N 



limit 18 



The leading (linear) corrections in a come from the three two-loop diagrams, shown in 
Fig. 4, which contribute to the polarization tensor. It turns out that all three contributions 
are antisymmetric in the indices \i and u, namely: 

TlW(p) . (25) 

This implies that there is no correction at this order to the longitudinal Fermion conductivty, 
&f,xx- The contribution to the Hall conductivty can be obtained by contracting each diagram 
in Fig. 4 with ^e' IuX p x and performing a trace over the Dirac spinor-components, which 



gives: 



tt(2) Q 4 6-D'-D f d °' k d ° 1 P \ ( k + P)Q \ (g + k ) (26] 

o{a+b) ° p2 V J (2n) D ' (2tt) d k 2 (k + p) 2 (k + q) 2 q 2 1 " 

n (2) = ? 9\ 6-d'-d f d °' k dD( l F(q,k,p) 

o(c) p 2 ^ J (2ir) D ' (2TT) D k 2 (k+p) 2 (k + q) 2 (k + q + p) 2 q 2 ' 1 ' 

F{q,k,p) = q 2 [3k 2 p 2 - 2{k ■ p) 2 + (k ■ p)p 2 - (k ■ p)(q ■ p) + (k ■ q)p 2 } 

-k 2 (q-p) 2 + (k ■ q)[2(k ■ p)(q ■ p) - 2(k ■ q)p 2 - (q ■ p)p 2 } . (28) 
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Since the integrals over k in Eqs. [27]) are convergent directly in D' = 3 dimensions, we 
perform them there to obtain 



n (2) _ £ U (3-D) [ d D q y _ p-(q + p) 
n o(a+b) 4p2 ^ J (2tt) d[ q\q + p\ q*\q + p\ 1 ' 



n (2) =^—u^- D ) f d q r q ' p ^ i p ' (l±A ) rsn N 

o(c) 4p 2^ J (2^) D[ q 2 {q + p) 2 q 2 \q + p\ 1 ' 1 ' 

The second terms in Eqs. Q29] , PO] ) are each logarithmically divergent in D = 3 but cancel 
one another exactly. The remaining finite contribution in D = 3 is finally: 



ni 2) W = ni 2 >(0) = n<$ a+b) + n g } = -^(i + ^) . (31) 



This result was obtained previously in Ref. |30| an d, up to an overall sign, independently 



in Ref. ^TJ (where the convention is that there is no imaginary unit in front of the Chern- 
Simons term). 

Upon inserting Eqn.(24) and ([31]) into (|13"D we obtain the final result for the Fermion 
conductivty at the anyonic Mott transition, valid up to first order in a: 

7T 

&F,xx = — , (32) 
o 

1 1 7T 2 

a F>xy = -- + -(l + —)a. (33) 

The contribution of -1/2 at a = comes from the second Dirac field which remains massive 
at the transition. The anyon conductivity follows by inverting to get the Fermion resistivity, 
and then using Eq. (14). Note again that at this order there is no correction to the 
longitudinal Fermion conductivity. An advantage of the 1/N method, which we describe in 
the next section, is that non-trivial corrections to <Jf,xx do appear at leading order in 1/N. 

To obtain the critical exponent i], let us consider the fermion self-energy S(p), defined 
as usual via the full Fermion propagator: 

S(p)- 1 = So^)- 1 - £(p) . (34) 
At one-loop order it is given by the diagram in Fig. 36: 
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£«(p) = -i|*. (35) 



The two- loop fermion self-energy diagrams are given in Fig. 5. In Ref. these diagrams 
have been evaluated to obtain: 

S (2) (P) = « ^^JZTd + + finite} . (36) 

Notice that at this order the renormalized mass remains zero, AM = E(p = 0) = 0, implying 
that conformal invariance survives quantum fluctuations at the transition point M = 0. The 
Fermion wave function renormalization constant is extracted to be 

x _ dS-\ P ) _ 1 

One can thereby obtain the anomalous dimension of the massless anyon field up to second 
order: 

^ ~ ~2 ~ "24^ ' { ) 

Finally let us consider the renormalization of the composite operator ipip, from which we 
can extract the critical exponent v. We insert this composite operator into one- and two-loop 
fermion self-energy diagrams, as shown in Figs. 4 and 5, with the insertion represented by a 
cross. To simplify the calculation, the external momentum associated with the ipip insertion 
is taken to be zero. We start with the one-loop diagrams, i.e. Fig. 6a, 66, and 6c. Within 
our regular ization scheme these diagrams are finite in D = 3: 

2 aiv\ A 

(6a) = (66) = -&r—?- , (39) 
o p 

(6c) = . (40) 

so that the anomalous dimension of the operator ipip indeed vanishes at one-loop order. 

However, the two-loop diagrams in Fig. 7 have logarithmic divergences. After lengthy 
but straightforward calculations we have (up to finite contributions which do not contribute 
to the anomalous dimension) 
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{7a 1 ) = 9 -^ D J d k 
(7a 2 ) + (7a 3 ) = ^ s - D J 
{7b x )=- 9 ^- D j 



(2tt) d k(k + p) 2 ' 
(2tt) d k(k + p) 2 ' 

g 4 3-d f d D k i 



(2tt) d /c(/c + p) 2 ' 

(7W + (76 3)^- D |^_T_. (41) 

The calculation of Figs. 7ci — 7c3 is a bit more complicated, because of the Dirac matrices. 
However, since we are only interested in the contribution proportional to the 2x2 unit 
matrix 1, it is sufficient to perform a trace over the product of Dirac matrices to obtain 
(again up to finite parts) 

-D'-D f d °'k d D q k 2 q 2 - (k ■ q) 2 



(7 Cl ) + (7c 2 ) + (7c 3 ) = 6sV 6 ~ D '- D J 



hV " i^n^- (42) 



(2tt) d ' (2tt) d k 2 q 2 ik + p) 2 {q + p) 2 {k + q + p) 2 
D q 1 
{2ti) d q 2 \q + p 

where we have used q ■ k = \[{q + k + p) 2 + p 2 — (q + p) 2 — (k + p) 2 ]. 
Finally, upon using the formula 

^V(SwW = i'3^ + <43) 
and putting everything together, we obtain the following divergent contribution to the 
Fermion self-energy with mass insertion: 

if °) = 1 + -^[-?- + InA] . (44) 
^ 16tt 2L 3-D p 2,i V ; 

Upon using the renormalization relation 

( r S^ 0) )r = Z i> Z W T< i£ } > ( 45 ) 

and Eqs. Q57D, (pB[), we can obtain the renormalization constant and the anomalous dimen- 
sion of the composite operator: 

y 1 1 (am 
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_ dZ^ _ 17g* 
Hi> - 9(1 j e } 24vr 2 ' y ' 



Finally, upon inserting Eqs. ( p8| ) and ([47]) into ([lq) , which defines the critical exponents r\ 
and v, we find up to second order in a: 

V = l~\a\ v-^l + ^a 2 . (48) 

3) 1/N Expansion 

In this subsection we consider a large N expansion which has the advantage of giving 
a non-vanishing leading correction to the Fermion longitudinal conductivity, in contrast to 
that found above. To this end, we consider generalizing the Lagrangian in Eqn. ([5]) for 
the massless Dirac field coupled to Chern-Simons gauge-field to include N massless Dirac 
fields. Once again we ignore for now the trivial contributions from the massive field. The 
appropriate Lagrangian is 

N g 1 

L = Wbvi 9 /* - i —nr T a v) + M ^ + h e ^X a ^ a X- ( 49 ) 
i=l vN 2. 

where as usual the coupling constant has been scaled by 1/N. The 1/N expansion is a formal 

summation of diagrams in powers oll/N, rather than powers of g 2 as in the previous section. 

Below we obtain the leading 1/N corrections to the conductivities at the transition and the 

critical exponents. 

At the critical point, M — 0, the Feynman rules in the Landau gauge now read 

So(p) = , (50) 
% p 

GS"(P) = , (51) 

rff(p) = i4gf ■ ( 52 ) 



However, in the 1/N expansion, rather than using the 'bare' gauge field propogator, Eq. 
(p)T|), it is more convenient to first sum up the chain of Fermion bubble diagrams shown in 
Figure. 6, since they each contribute at the same order, namely O(l/N ). (Each fermion 
loop now carries an extra factor N because of the summation over flavors, and the two 



interaction vertices carry g 2 /N.) Upon using Eqs. ([12]) and ([24]) to help sum up the bubbles 
in Fig. 8, we obtain an effective (or dressed) gauge-field propagator, denoted by a cross line 
in the Figures, as follows 

GT ff (p) = A d P / P +B € —^, (53) 

A = ^ , B = J-g— . (54) 

A feature of the 1/N expansion in the present model is that non-trivial corrections to the 
critical exponents arise already at leading order (in 1/N), in contrast to the weak coupling 
expansion, where it was necessary to go to two-loop order to see corrections. In addition, as 
verified by explicit calculations below, the 1/N expansion maintains conformal invariance 
at the transition point M = 0, at least at the leading order in 1/N. 

To leading order, 0(1/N), the Fermion self-energy is given by the diagram in Fig. 9. 
We have verified that it does not shift (or renormalize) the Fermion mass from M = 0. To 
extract the term that is proportional to we perform $Tr i> on the diagram in Fig. 9. 
Keeping only the divergent part, we thereby find 

^"'^ih +ln 4 )] - (55) 

which gives directly the anomalous dimension: 

^ = -~6^N ■ (56) 
Upon insertion of the operator ipif) into the one-loop Fermion self-energy, we have the 
three diagrams in Fig. 10, which up to finite parts give a contribution: 

d D k 1 



(l0a) = 2A 9 -^ D J 



(2tt) d (k + p) 2 k 

9 A / A 2 r>2\ . 3-D f d D k 1 



m + (!0c) = -foA> - B V -o j JLJL^ . (57) 



Performing the loop integrals then gives 



(^-BV , 9 2 
L 3^ 2 J 27r 2 iV3-D 



^ = 1 - ~ " TT^^r— , (59) 
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and the anomalous dimension 



$ (A 2 -B 2 )g\ 2 



Finally, upon inserting Equations flog ) and (|60|) into (pf) we obtain the leading 1/N 
expressions for the critical exponents: 

, 8 a 2 1 

V = 1 " 3 64TJ^N ' (61) 
128 (128 -(7ra) 2 ) a 2 1 

3 (64 + (vra) 2 ) 2 iV ' 1 J 

Next we calculate the order 1/iV corrections to the conductivities at the critical point. 
The relevant Feynman diagrams which contribute to the Fermion polarization tensor are 
listed in Figs. 11 and 12. The diagram in Fig. 12 is zero by Furry's theorem, because it 
contains a closed Fermion loop attached to an odd number of gauge-field lines. (Essentially 
this is a consequence of charge-conjugation invariance, since the Chern-Simons gauge boson 
has odd charge parity.) Fig. 11a, 116, and 11c are linear in and therefore linear in e fluX 
and in (5 tJ,u — p^p u /p 2 ). Thus the odd contribution, U (p), can be read off directly from (|3T|): 

9 4 B , vr 2 1 
1 (1 H )_ 

16tt 21 4 >N 



^o ia+b+c) (p) = ^- 2 (l + -)- T - (63) 



Calcualting the even contribution, U e (p), however, is much more complicated. In addition to 
the cumbersome Dirac trace, the evaluation of the Feynman integrals are highly non-trivial. 
We leave the details to an appendix and here only quote the final result, 

It should be emphasized that all divergent contributions to the polarization tensor 11^ (p) 
cancel, so that in the 1/N expansion the /3-function /3(g) vanishes, and conformal invariance 
survives quantum fluctuations at the transition point. 

Finally we can obtain the Fermion conductivities at the critical point from the above 
polarization tensors. Eq. (13) can be used to obtain the conductivity per flavor of massless 
Dirac field. To obtain the total Fermion conductivity per flavor we must add to (yp xy the 
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contribution of -1/2 from the massive Dirac field. The final result to leading order in 1/N 
is: 

= -\ + led + t 'mTT^FF • (65) 

7r 12a 2 1 . 

8 64 + {na) z TV 

IV. DISCUSSION 

As we have seen, the critical properties of the Mott anyon transition vary continuously 
with the anyon statistics. Thus the model describes a line of fixed points which are char- 
acterized by the statistics parameter a, which varies from Fermions, a = 0, to Bosons, 
a — 1. In Section III we calculated the critical exponents and universal conductivities at 
the anyon Mott transition as an expansion around the Fermion point. Specifically, in terms 
of the deviation from Fermi statistics, a, we obtained critical exponents up to second order 
and conductivities to first order. The critical properties at the Mott transition can also be 
obtained at a = 1 directly in terms of a Bosonic scalar (U(l)) field theory, or equivalently 
the 3D XY model fllfl . The exponents are of course known quite accurately for the 3D 
XY model, and recently an estimate for the conductivity has been obtained from Monte 
Carlo simulations and a large N expansion ||12||. Also, within a large N calculation Wen and 



Wu |13[ have recently performed an expansion for a near one by coupling a Chern-Simons 
gauge-field to the U(l) Boson field. It is instructive to compare these various results for 
the critical properties of the Mott transition in order to see the trends in exponents and 
conductivities as one varies the particle statistics from Fermi to Bose. As we discuss below, 
the observed trend can give one some insight into possible connections between transitions in 



the integer quantum Hall effect |35| , and the magnetic field-tuned superconductor insulator 
transition in thin films |36 |, [37]. 



Consider first the exponents rj and v. In Section III we found that to second order in a, 

77(a) = 1- -a 2 + 0(a 3 ) , v-\a) = 1 + — a 2 + 0(a 3 ) . (67) 
3 6 

For the 3D XY model, which corresponds to a = 1, the exponents are given by roughly, 
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7j(a = 1) ~ 1/50, u~ 1 (a = 1) ~ 3/2 . (68) 

Notice that both 77 and u are smaller for Bosons than for Fermions. Moreover, the expansion 
from the Fermion end indicates that the initial deviations for small a are to reduce the 
exponents, suggesting that 77(a) and v(pt) might be monotonically decreasing functions of 
a. Unfortunately, as is clear from this expansion, the leading order term does not give a 
reliable estimate at a = 1. This should be contrasted with the e = 4 — D expansion for the 
XY model @ which gives reasonable exponent values for the 3D model when the low order 
results are extrapolated to e = 1. 

It is also instructive to compare the 3D XY exponents with the exponents obtained from 
the large N Fermion expansion in Section 3: 

8 a 2 1 

^ = 1 - 364 +(7 ra) 2 iV' (69) 
128 (128 -(7ra) 2 ) Q 2 1 

3 (64 + (Tra) 2 ) 2 N ' { ' 

If we put a = 1 these become 

v (a = 1) = 1 - (0.036.. .)(l/iV) + 0(l/iV 2 ) , (71) 
v-^a = 1) = 1 + (0.924...) (1/iV) + 0(l/iV 2 ) , (72) 

which, when extrapolated to N=l, should become equal to the 3D XY model exponents. 
Once again, although the leading term in the expansion has the 'correct' sign, the extrapo- 
lations to N=l using only the first term are clearly not very reliable. 

Next we consider the universal conductivities at the Mott transition. For Fermions, 
a = 0, we obtained in Section III that a xx = tt/8 and a xy = —1/2. For the Boson Mott 



transition (in zero field) a xy = and the recent numerical estimates |12j give a xx (a = 1) = 
0.285 ± 0.02, a value slightly smaller than in the Fermion case. It is instructive to see if the 
expansion about the Fermion point for small a gives the correct trends. For this purpose it 
is both more natural and easier to compare resitivities, since as (14) shows the longitudinal 
anyon resistivity is simply equal to the Fermion longitudinal resistivity. Upon inverting the 
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perturbative results in Equation (32) and (33) for the Fermion conductivity, we obtain, using 
(14), an expansion for the anyon longitudinal resistivity: 

Pxx(a) = (0.9715...) + (2.084...)a + 0(a 2 ) . (73) 
Simarlarly, the transverse anyon resistivity is given by: 

p xy (a) = (1.237...) - (0.492.. .)a + 0(a 2 ) . (74) 
The universal resistivities at a = 1 follow from the Boson Monte Carlo simulations and are 







p xx (a = 1) = 3.5 ±0.2, Pxy ( a = 1) = 0. (75) 

Notice that the sign of the leading small a correction to both p xx and p xy are such that 
the values tend towards the boson values in (75). This result suggests that, just as with 
the critical exponents, the components of the universal resistivity tensor vary monotonically 
upon moving along the fixed line of critical points from Fermion to Boson statistics. 

One might be tempted to compare directly our large-N results to those obtained in Ref. 
[ 13]] for N-flavors of Bosons coupled to the Chern-Simons field. However, the two large 



N-theories are probably not continuously connected to one another (upon varying a) since 
the extension to N fields brings into the Lagrangian of the theory a new SU(N) symmetry, 
which is spontaneously broken upon crossing the Bosonic transition, while unaffected across 
the Fermionic phase transition. 

We now turn to the relevance of the results obtained in this paper to two-dimensional 
experimental systems which exhibit zero temperature quantum phase transitions. Unfortu- 
nately we cannot make direct contact with experiment, since our simplified lattice model 
ignores both disorder and long-ranged Coulomb interactions. Nevertheless, the notion of 
a fixed line of critical points parameterized by the statistics of the underlying particles, 
along which the exponents vary continuosly, and monotonically, is presumably rather more 
general. 
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To be specific, consider a more realistic model of anyons, with a long-ranged Coulomb 
interaction moving in a quenched random potential and external magnetic field. In the 
Boson limit, a = 1, as parameters are varied this model presumably undergoes a transition 
from a superconducting phase to a localized Bose glass insulator. This transition is believed 
to be in the appropriate universality class for real disordered superconducting films, which 
are tuned with an external magnetic field from the superconducting into insulating phases 
]37fl . This magnetic field tuned superconductor- insulator transition has been recently studied 
experimentally in considerable detail []36| |. In the Fermion limit, a = 0, with strong disorder 



the model is presumed to exhibit transitions between integer quantum Hall plateaus, for 
which there is also considerable experimental data |35] . The critical properties of these two 



experimentally accessible phase transitions are thus presumably end points of a fixed line of 
critical points which interpolates between them. Moreover, transitions between fractional 
plateaus in the quantum Hall effect are probably described correctly by this same model 
with fractional statistics a. For example, the transition between the so-called Hall insulator 
and the v = 1/3 quantum Hall plateau can be described as a condensation of fractional 
statistics particles with a = 2/3. 

Although at present we cannot calculate analytically (or numerically) the critical prop- 
erties along the fixed line of these disordered anyon models, it is instructive to compare the 
critical properties measured experimentally with the trends obtained in our simplified clean 
lattice model. For example, experiments [|35| on the transition between integer plateaus in 



the Hall effect find that the correlation length exponent, u, times the dynamical exponent, 
z, to be given by roughly 7/3. Since one expects z=l in the presence of 1/r Coulomb inter- 



actions this gives an estimate for v ~ 7/3. At the magnetic field tuned superconductor 
insulator transition, on the other hand, with z=l assumed, the experiments give a much 
smaller value, v ~ 5/4. Thus upon moving along the presumed line of fixed points from the 
Fermion to Boson end, the exponent v apparently decreases. It is noteworthy that this same 
trend, a decreasing of v moving from Fermion to Boson, was what we found along the fixed 
line of Mott anyon transitions. Perhaps it is generic that Fermion transitions are closer to 
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their lower critical dimension, with a larger u, than their Boson counterparts. 

It is also amusing to compare the experimental values for the universal resistivites at the 
two transitions. At the field tuned superconductor-insulator transition [36[] the longitudinal 
resistivity is found to cluster in the range between 0.8 and 1.0, in units of h/(2e) 2 , where the 
Cooper pair has charge 2e. At the transition between integer plateaus in the quantum Hall 
effect, experiments find |39|] values of a xx of roughly 0.2. Combining this with o xy = 1/2 gives 



p xx of roughly 0.7, a value slightly smaller than at the superconductor-insulator transition. 
It is interesting that at the anyon Mott transition we also find the longitudinal resistivity 
increases moving from the Fermion to Boson end. It is of course unclear whether or not this 
trend is a generic property of anyon phase transitions. 

Finally, it is worth re-empasizing that the existence of a fixed line of critical points 
characterized by the anyon statistics in our simple lattice anyon model, underscores the 
importance of the fluctuating gauge field. At a critical point, or more generally in a gapless 
phase, the fluctuating Chern-Simons gauge field cannot be simply thrown away, as in the 
'flux-smearing' mean field approaches. More specifically, we expect that the presence of this 
fluctuating gauge field will most likely make the transition between plateaus in the fractional 
quantum Hall effect in a different universality class from that between integer plateaus. 
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APPENDIX: Calculation of U e (p) at order 0(1/N) 

In this appendix we present the details for the two-loop calculation of the even part of 
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the gauge-boson self-energy Tl e (p) in the 1/N expansion. 
Some useful formula for the 2x2 7-matrices are 



-3 - 1 , 



a A u A 

7 7 7 =7, 



7^7 V = 25 aX l - TT , 
r r ^ J VYe^ T p T = 2p n l. 



(76) 
(77) 
(78) 
(79) 
(80) 
(81) 
(82) 



The relevant diagrams are those in Fig. 11. Note that U e = ±5^ U W V '. We have Il'f+V 
and IT^ are therefore given by 



A£l_ 6-D'-D 



dP'k d D q 7^+^)7^ #7"(#+ ^)7 CT #(<5 CT V - gV) 



(2tt) b ' (2tt) 



1? 



Tr 



k 4 (k+p) 2 (k + g) 2 g 3 



.Ag A q_ d >_ d f d D 'k d D q q ■ (q + k)(2k ■ pk ■ q + k 2 k ■ q — k 2 p ■ q) 



p N 



and 



(2tt) d ' (2tt) d 



k\k + p) 2 (k + q) 2 q 3 



(83) 



A 9 6-D'-D 

' p2N^ 

Ag r-d'-d 

p JM 



d D 'k d D q 7^+ 4+ 4)l u ft(5 av q 2 - q a q v ) 



(2tt) d ' (2n) D 
d D 'k d D q 



Tr 



k 2 (k + p) 2 (k + q) 2 (k + q + p) 2 q 3 
G(k,q,p) 



(2tt) d ' (2tt) d k 2 (k + pf{k + q) 2 (k + q + p) 2 q 3 



(84) 



respectively, where 



G(k, q,p) = k 2 [q A — q 2 p ■ q — 2q 2 k ■ q — 2q 2 k ■ p — (g ■ p) 2 + q 2 p 2 — k 2 q 21 

+ (k ■ q)[ — 2q 2 k ■ q + 2(q ■ p){k ■ p) — 2q 2 k ■ p — 2q 2 q ■ p — p 2 k ■ q] 
+q 2 (k ■ p)\q 2 — q ■ p — 2k ■ p] . 



(85) 



Performing the convergent ^-integrations in Eqs. (83, |84|) at D' = 3, we have 
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A g A 3 _ D f d D q r 2g ■ p + p 2 q-p 



TT g ,,3--D / I" g ^ 1 /Q«\ 

ll e (o+6) = — ^rr/i / / \ n [ 3] — ; — i — ; — r • ^ b > 

p 8l\ J {1-k) u q A \q + p\ q \q + p\p 



^ A g 4 3 _ D r d D q 2q ■ p + p 2 Ap 1 2g 4 + 5p 2 q 2 + 2p 4 

e ^ p 8N^ J (2tt) d q 3 \q + p\ Q 2 \q + p\ q\q + p\ (q ' P) ( i 2 \ < i + p\p 



Due to the factor q-p in the denominator of the last term in (57), it is hard to obtain an 
analytic expression with the dimensional regularization. So we will set D = 3 and introduce 
a momentum cut-off A in both (|5BJ) and (|57D. We will see the logarithmically divergent 
terms indeed cancel. Also we will discard linearly divergent terms, whose appearance is an 
artifact of this regularization by momentum-cutoff. (Such divergences do not appear in a 
regularization that preserves Lorentz invariance.) This results in a finite contribution to IT e 
at the order 0(1/N): 

Ag 4 f A d 3 q q-p Ap 1 2g 4 + 5p 2 q 2 + 2p 4 

ll e ( a+b+c ) — / — ; — + — ; 1 + ; 1 + — ; — ; ; . (88 ) 

p oJM J (27rj d q \q + p\p q \q + p\ Q\Q + P\ {q ■ p)q \q + p\p 

Except for a factor — ^7^x2, the first three terms are, respectively, 

2 1 A 2 A 2 

In—, 8 + Aln— 

9 3 p 2 p 

The fourth term is a bit messy: 



In— , 8 + Aln— , - 2 . (89) 



1 yi 2 + by + 2y 2 1 VT+y^/y 1 r^r 2 + 5y + 2y 2 ; v / T+y-l 



2 7 yv/TTy VT+^+vf 2 7i t/v^ Vl+y+l 
1 /■ 1 . 2 + 5y + 2y 2 VT+y-y^ c 22 A 2 

2io ^ y^T+y /n yr+^ + ^ + C 6% 2 ' (9 ° } 



where 



<7 = -51 - l -[ln4±lf + lM Zn ^±l = 2 . 137820 914 . (91) 
9 2 L v/2-l J 3 V2-1 

The integral in ( |90f) is numerically evaluated to be —6.405956452 • • • . The logarithmic 
divergences cancel out and the final finite result is 

A o 4 3A o 4 

D.M = -1.509642238-^ w ___JL_ . ( 92) 



27 



REFERENCES 

[1] F. Wilczek, Phys. Rev. Lett. 48, 1144 (1982); Phys. Rev. Lett. 49 957 (1982); Y.S. Wu, 
Phys. Rev. Lett. 53, 111 (1984). 

[2] A. Fetter, C. Hanna and R. Laughlin,Phys. Rev. B39 (1989) 9679. 

[3] S.M. Girvin and A.H. Macdonald, Phys. Rev. Lett. 58 1252 (1987). 

[4] S.C. Zhang, T.H. Hansen and S. Kivelson, Phys. Rev. Lett. 62, 82 (1989); D.H. Lee 
and S.C. Zhang, Phys. Rev. Lett. 66, 1220 (1991). 

[5] N. Read, Phys. Rev. Lett. 62, 86 (1989). 

[6] J.K. Jain, Phys. Rev. Lett. 63 (1989) 199; Phys. Rev. B40 (1990) 8079. 

[7] B. Blok and X.G. Wen, Phys. Rev. B42 (1990) 8133; X.G. Wen and A. Zee, Nucl. 
Phys. B351 (1990) 135; X.G. Wen, Mod. Phys. Lett. B5 (1991) 39. 

[8] B.I. Halperin, P.A. Lee and N. Read, preprint, 1992. 

[9] D.H. Lee, S. Kivelson and S.C. Zhang, Phys. Rev. Lett. 68, 2389 (1992); preprint, 1992. 

[10] G.S. Canright, S.M. Girvin and A. Brass, Phys. Rev. Lett. 63, 2295 (1989). 

[11] M.P.A. Fisher, P.B. Weichman, G. Grinstein and D.S. Fisher, Phys. Rev B40, 546 
(1989). 

[12] M. Cha, M.P.A. Fisher, S.M. Girvin, M. Wallin and A.P Young, Phys. Rev B44, 6883 
(1991). 

[13] X.G. Wen and Y.-S. Wu, MIT preprint (July, 1992). 

[14] A tight binding model with half of a flux quantum per plaquette arises in 'flux-phase' 
mean field treatments of frustrated two-dimensional antiferromagnets. See, for example, 
X.G. Wen et. al. Phys Rev. B39, 11413 (1989). 



28 



[15] D.J. Thouless, M. Kohmoto, N.P. Nightingale and M. den Nijs, Phys. Rev. Lett. 49, 
405 (1982). 

[16] We are grateful to R. Shankar for the following particularly succinct derivation of the 
final Dirac representation in Equation (7). 

[17] W. Siegel. Nucl. Phys. B156 (1979) 135; J.J. Schonfeld, Nucl.Phys. B185 (1981) 157; 
R. Jackiw and S. Templeton, Phys. Rev. D32 (1981) 2291. 

[18] M.P.A. Fisher, G. Grinstein and S.M. Girvin, Phys. Rev. Lett. 64, 587 (1990); X.G. 
Wen and A. Zee, Int. J. Mod. Phys.B 4, 437 (1990). 

[19] D. Amit, Field Theory, Renormalization Group and Critical Phenomena, (Second Edi- 
tion), World Scientific, 1988. 

[20] G. Parisi, Statistical Field Theory, Addison- Wesley, 1988. 

[21] P. Ramond, A Priemer of Field Theory, Addison- Wesley, 1988. 

[22] S. Deser, R. Jackiw and S. Templeton, Phys. Rev. Lett. 48 (1982) 975; Ann. Phys. 
(N.Y.) 140 (1982) 372. 

[23] S. Coleman and B. Hill, Phys. Lett. 159B (1985) 184. 

[24] N. Redlich, Phys. Rev. D29 (1984) 2366. 

[25] A. Niemi, GW. Semenoff and Y.-S. Wu, Nucl. Phys. B276 (1986) 173. 

[26] GW. Semenoff, P. Sodano and Y.-S. Wu, Phys. Rev. Lett. 62 (1989) 715. 

[27] W. Chen, G.W. Semenoff and Y.-S. Wu, Mod. Phys. Lett. A5 (1990) 1833; in Physics, 
Geometry and Topology, Proc. of Banff Summer School on Particles and Fields, August 
1989; (Plenum Pub. Cor.), p.553, 1990. 

[28] W. Chen, G.W. Semenoff and Y.-S. Wu, Phys. Rev. D44 (1991) R1625; and Phys. Rev. 
D, in press. 

29 



[29] W. Chen, Phys. Lett. 251B (1990) 415. 
[30] W. Chen (unpublished) (1990). 

[31] V.P. Spiridonov and F.V. Tkachov, Phys. Lett. B260 (1991) 109. 

[32] There is a huge literature for perturbative calculation of the induced Chern-Simons term 
in various D = 2 + 1 quantum field theories other than Abelian Chern-Simons coupled 
to non-self-interacting matter. An incomplete list is: Y. Kao and M. Suzuki, Phys. rev. 
D31 (1985) 2137; M. Bernstien and T.J. Lee, Phys. Rev. D32 (1985) 1020; R. Pisarski 
and S. Rao, Phys. Rev. D32 (1985) 2081; T. Lee, Phys. Lett. 171B (1986) 247; V.P. 
Spidonov, JETP Lett. 52 (1990) 1112; L.V. Avdeev, G.V. Grigoriev and D.I. Kazakov, 
CERN preprint TH-6091/91 (1991); G. Ferretti and S.G. Rajeev, Mod. Phys. Lett, bf 
A7 (1992) 2087; S.H. Park, Phys. Rev. D45 (1992) R3332; D.K. Hong, T. Lee and 
S.H. Park, Korea preprint SNUTP 92/91 (1992); J. Chay, D.K. Hong, T. Lee and S.H. 
Park, Korea preprint SNUTP 92/92 (1992). 

[33] In fact the beta-function /3(g) vanishes by explicit calculation at two loops even in non- 
Abelian theory [p^j ; it is proved to vanish up to all orders in perturbation theory in 
massive Dirac theory, and is conjectured to do so in massless theory. An argument for 
the latter case can be given as follows: The Chern-Simons Lagrangian is gauge invariant 
only up to total divergence terms. But on the other hand, only strictly gauge invariant 
local counterterms can arise in a gauge-invariant regularization and renormalization 
scheme. So no infinite counterterm for the Chern-Simons Lagrangian can appear. 

[34] The divergence of the pole 2/e = 2/(3 — D) in the physical dimension D = 3 corresponds 
to ln(K 2 / /i 2 ) in the regularization with momentum cutoff A. 

[35] H.P. Wei, D.C. Tsui, M.A. Paalanen and A.M.M. Pruisken, Phys. Rev. Lett. 61, 1294 
(1988). 

[36] A.F. Hebard and M.A. Paalanen, Phys. Rev. Lett. 65, 927 (1990). 



30 



[37] M.P.A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 
[38] M.E. Fisher, Rev. Mod. Phys. 46, 597 (1974). 
[39] D.C. Tsui (private communication). 



31 



FIGURES 




Fig. 1 The four sites in the unit cell of the two-dimensional square lattice tight binding 
model with one-half of a flux quanta per plaquette and a staggered periodic potential. 
Near neighbors have hopping strength t except the bold line which is —t. The next- 
near neighbors have a hopping strength it'/ 4, and the staggered on-site potential has 
strength m. 




Fig. 2 The full current- current correlation for the external electromagnetic field. The 
black spots stand for the currents. The striped disk represents the exact Chern-Simons 
self-energy and the dashed line the Chern-Simons propagator. 
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a b 

Fig. 3 (a) Chern-Simons and (b) fermion self-energies at 0(e 2 ). The solid line stands 
for the fermion propagator. 
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Fig. 5 Fermion self-energy at 0(e 4 ). 





Fig. 6 Vertices with composite operator insertions at C(e 2 ). The cross stands for the 
composite operator tptp. 
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^3 b 3 C 3 

Fig. 7 Fermion vertex with the insertion of the composite operator ijjip at 0(e 4 ). 




Fig. 8 Summation of infinite series of one-loop Fermion bubble chains gives the dressed 
Chern-Simons propagator, which is of order O(N ). The cross line is the dressed 
Chern-Simons propagator, the solid lines are the Fermion propagators, the dashed 
lines are the bare Chern-Simons propagators. 
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Fig. 9 Fermion self-energy at 0(1/N). 



+ + + 7\ + 

+ + + + 

a b c 

Fig. 10 Fermion vertex with an insertion of the composite operator at 0(1/ N). 
The cross stands for the composite operator 'ipiipi- 
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a b 
Fig. 11 Chern-Simon self-energy at 0(1/ N). 

<x> 

Fig. 12 Null diagram at order 0(1 /N). 



37 



